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Abstract 

We present a detailed analysis of the neutrino-induced muon signals coming from 
neutralino pair-annihilations inside the Sun and the Earth with particular empha- 
sis for light neutralinos. The theoretical model considered is an effective MSSM 
without gaugino-mass unification, which allows neutralinos of light masses (below 50 
GeV). The muon events are divided in through-going and stopping muons, using the 
geometry of the Super-Kamiokande detector. In the evaluation of the signals, we 
take into account the relevant hadronic and astrophysics uncertainties and include 
neutrino oscillation and propagation properties in a consistent way. We derive the 
ranges of neutralino masses which could be explored at neutrino telescopes with a 
low muon-energy threshold (around 1 GeV) depending on the category of events and 
on the values of the various astrophysics and particle-physics parameters. A final 
analysis is focussed to the upward muon fluxes which could be generated by those 
neutralino configurations which are able to explain the annual modulation data of the 
DAMA/LIBRA experiment. We show how combining these data with measurements 
at neutrino telescopes could help in pinning down the features of the DM particle and 
in restraining the ranges of the many quantities (of astrophysics and particle-physics 
origins) which enter in the evaluations and still suffer from large uncertainties. 
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1 Introduction 



In the papers of Ref. [1] it was shown that hght neutrahnos with a mass in the range 7 GeV 
< < 50 GeV are interesting candidates for particle dark matter, with direct detection 
rates accessible to experiments with large exposure (of order 100 000 kg day) and low 
energy threshold (of a few keV). This population of light neutralinos arises in the Minimal 
Supersymmetric extension of the Standard Model (MSSM) when the unification of gaugino 
masses at a Grand Unified (GUT) scale is not assumed [2]. Indeed, in this supersymmetric 
framework the lower bound of ~ 7 GeV on the neutralino mass is set by a cosmological 
bound on the neutralino relic density [1]. This is at variance with the lower bound > 
50 GeV, which is derived from the LEP lower limit on the chargino mass within the MSSM 
with gaugino mass unification at the GUT scale. 

A direct comparison of the theoretical predictions of Ref. [1] with experimental data was 
made possible when the DAMA Collaboration published the results of its measurements 
collected with a Nal detector of 100 kg over 7 annual cycles [3]. Actually, in Ref. [1] it was 
proved that the population of light neutralinos fitted well these DAMA results. 

The inclusion of the channeling effect in the experimental analysis [5] and the further 
results of the DAMA/LIBRA combined data [6] confirmed the good agreement of the 
predictions of Ref. [1] with the experimental data, as discussed in Refs. [7, 8]. 

Possible indirect effects of light neutralinos were studied in Refs. [8, 9], where it was 
shown that present cosmic antiproton data set sizable constraints on the supersymmetric 
space; measurements of cosmic antideuterons [10] with forthcoming airborne experiments 
[11] were indicated in [8, 9] as a very promising investigation mean for the light neutralino 
population. 

In Ref. [9] also possible signals induced in neutrino telescopes by the neutrinos produced 
by pair-annihilations of light neutralinos captured in the Earth and the Sun were discussed. 
The analysis was carried out in the case of a single category of events, those due to upward 
through-going muons. 

In the present paper we resume the analysis of Ref. [9], by implementing and extending 
it in various distinctive features: a) the particle-physics uncertainties in hadronic quan- 
tities which affect the capture rate of relic neutralinos by the celestial bodies are taken 
into account and discussed in detail; b) all main processes that occur during the neutrino 
propagation: neutrino oscillations and neutrino incoherent interactions with matter are 
included; c) the investigation comprises also the category of upward stopping muons (ac- 
tually, this turns out to be the most promising case for upward-going muons induced by 
pair annihilation of light neutralinos). 

Indirect evidence of WIMPs in our halo by measurements of upward-going muons at 
neutrino telescopes, generated by WIMP pair-annihilation in celestial bodies, has been the 
subject of many investigations in the past, see for instance [12, 1.3, 14, 15]. 

Recently, a number of papers appeared where possible signals at neutrino telescopes 
due to pair annihilation of light DM particles are discussed. These consider either generic 
WIMPs with assumed dominance of specific annihilation channels [16, 17] or discuss specific 
DM candidates other than neutralinos: WIMPless dark matter and mirror dark matter 
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[17], leptonically interacting dark matter [IS] and DM particles which directly annihilate 
to neutrinos [19]. All these investigations are restricted to signals from the Sun. 

At variance with these analyses, the present paper deals with neutrino signals produced 
in cosmic macroscopic bodies by pair annihilation of light neutralinos whose specific super- 
symmetric properties are those which satisfy all existing particle-physics and cosmological 
constraints, as analyzed in Refs. [1, <s]. Furthermore, our investigation entails neutrino 
signals expected from the Sun as well as those from the Earth. 

Our results are given first for the whole population of light neutralinos, then for those 
subsets which are selected by the DAMA/LIBRA results [C], when these are interpreted 
in terms of relic neutralinos. We separate the case where the channeling effect is included 
from the one where this effect is neglected. 

The scheme of the present paper is the following. In Sect. 2 the main features of the 
super symmetric model employed here are summarized. The formulae providing the capture 
rates of relic neutralinos by celestial bodies and the annihilation rates due to neutralino 
pair-annihilation are recalled in Sect. 3, together with some properties of the hadronic 
quantities which enter in the evaluations of the neutralino-nucleon cross section. The 
generation and the propagation of the neutrino fluxes is described in Sect. 4, their ensuing 
muon fluxes are derived in Sect. 5. Results and conclusions are given in Sect. 6. 

2 Theoretical model 

The supersymmetric scheme we employ in the present paper is the one described in Refs. [1, 
8]: an effective MSSM scheme (effMSSM) at the electroweak scale, with the following 
independent parameters: Mi, M2, yU, tan/3, m^, 'mg, m^- and A. Notations are as follows: 
Ml and M2 are the U(l) and SU(2) gaugino masses (these parameters are taken here to be 
positive), /i is the Higgs mixing mass parameter, tan/3 the ratio of the two Higgs v.e.v.'s, 
niA the mass of the CP-odd neutral Higgs boson, rriq is a squark soft-mass common to all 
squarks, mj is a slepton soft-mass common to all sleptons, and A is a common dimensionless 
trilinear parameter for the third family, = A^ = Arrig and = Arrij (the trilinear 
parameters for the other families being set equal to zero). In our model, no gaugino-mass 
unification at a Grand Unified (GUT) scale is assumed. The lightest neutralino is required 
to be the lightest supersymmetric particle and stable (because of R-parity conservation). 

The numerical analysis presented in the present paper was performed by a scanning of 
the supersymmetric parameter space, with the following ranges of the MSSM parameters: 
1 < tan/3 < 50, 100 GeV < < 1000GeV,5GeV < Mi < 500 GeV, 100 GeV < M2 < 
1000 GeV, 100 GeV < nig, mj < 3000 GeV, 90 GeV < < 1000 GeV, -3 < A < 3. 

The supersymmetric parameter space has been subjected to all available constraints 
due to accelerator data on supersymmetric and Higgs boson searches (CERN e~^e~ collider 
LEP2 [20] and Collider Detectors DO and CDF at Fermilab [21]) and to other particle- 
physics precision results: measurements of the 6 — > s -|- 7 decay process [22] [2.89 < 
B{b — > s -|- 7) ■ 10~^ < 4.21 is employed here], the upper bound on the branching ratio 
5i?(5° ^ /i~ + /i+) [23] [we take BR{B'^ ^ /i" + /i+) < 1.2 ■ 10"^] and measurements 
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of the muon anomalous magnetic moment = {g^ — 2)/2: for the deviation Aa^ of the 
experimental world average from the theoretical evaluation within the Standard Model we 
use here the range -98 < Aa^ ■ 10^^ < 565. For other details concerning these constraints 
see Refs. [1, 8]. 

Also included is the cosmological constraint that the neutralino relic abundance does 
not exceed the maximal allowed value for cold dark matter, i.e. VL-^h? < {0,cDMh'^)ma.x- 
We set (ficDM^^)max = 0.122, as derived at the 2a level from the results of Ref. [24]. We 
recall that this cosmological upper bound implies on the neutralino mass the lower limit 
'^x ^ 7 GeV [1]. 



3 Capture rates and annihilation rates 

In the present section we summarize the main formulae which we employed to evaluate the 
capture rates of neutralinos by celestial bodies and their pair-annihilation in neutrinos. 

The capture rate C of the relic neutralinos by a macroscopic body is given by the 
standard formula [25] 

where is the neutralino mean velocity, a^^i is the cross section of the neutralino elastic 
scattering off the nucleus i of mass m,, Msfi is the total mass of the element i in the 
body of mass Mb, {Vesc)i is the square escape velocity averaged over the distribution of 
the element i, Xi is a factor which takes account of kinematical properties occurring in the 
neutralino-nucleus interactions. 

The neutralino local density is denoted by p^. One can assume that is equal to the 
local value of the total DM density po, when the neutralino relic abundance {Q^h^) turns 
out to be at the level of a minimal (f^cDAf ^^)min consistent with po- On the contrary, when 
is smaller than {^cDMh'^)mm, the value to be assigned to has to be appropriately 
reduced. Thus we evaluate fl^h'^ and we determine p^ by adopting a standard rescaling 
procedure [12]: 



Px = Po , when Q.^h'^ > {VLcDuh^^ 



/mm 



VL (2) 

Px = POTF> ^~r2\ ' ^X^^ < (^CDA/^^)mm 

\}l'CDMil' jmin 

Here (i^cDM^^)min is set to the value 0.098, as derived at 2a level from the results of 
Ref. [24]. The neutralino relic abundance is evaluated with the procedure discussed in 
Ref. [1]. 

As for the velocity distribution of relic neutralinos in the galactic halo we use here, 
for definiteness, the standard isothermal distribution parametrized in terms of the local 
rotational velocity vq. It is however to be recalled that the actual distribution function 
could deviate sizably from the isothermal one (see, for instance, Ref. [26] for a systematic 
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analysis of different categories of distribution functions) or even depend on non-thermalized 
effects [27]. Also the possible presence of a thick disk of dark matter could play a relevant 
role in the capture of dark matter by celestial bodies [28]. The local rotational velocity 
Vq is set at three different representative values: the central value Vq = 220 Km s~^ and 
two extreme values Vq = 170 Km s~^ and vq = 270 Km s~^ which bracket the Vq physical 
range. Associated to each value of vq we take a value of po within its physical range 
established according to the procedure described in Ref . [26] . In conclusion, we will provide 
the numerical results of our analysis for the following three sets of astrophysical parameters: 
1) Vq = 170 Km s-\ Po = 0.20 GeV cm'^; 2) vq = 220 Km s-\ po = 0.34 GeV cm'^; 3) 
Vq = 270 Km s~^, Po = 0.62 GeV cm~^. Note that these values of po correspond to the 
case of maximal amount of non halo components to DM in the galaxy [26]. 

The annihilation rate Fann is expressed in terms of the capture rate by the formula 



Tann = ^tanh^ ^-^ j , (3) 

where t is the age of the macroscopic body (t = 4.5 Gyr for Sun and Earth), ta = 
{CCa)~^^'^, and Ca is the annihilation rate per effective volume of the body, given by 

_ <av> l' \3/2 

In the above expression Vq is defined as Vq = {3mpiT/{2p x 10 GeV))^/^ where T and p 
are the central temperature and the central density of the celestial body. For the Earth 
(T = 6000 K, p = 13 g ■ cm-3) Vq = 2.3 x lO^^cm^, for the Sun (T = 1.4 x 10^ K, 
p = 150 g ■ cm~^) Vo = 6.6 x 10^^ cni'^. cr is the neutralino-neutralino annihilation cross 
section and v is the relative velocity. The thermal average < cri; > is calculated with all 
the contributions at the tree level as in Ref. [30], with the further inclusion here of the two 
gluon annihilation final state [13]. 

We recall that, according to Eq. (3), in a given macroscopic body the equilibrium bet- 
ween capture and annihilation {i.e. Fann ~ C/2 ) is established only when t > ta- It is 
worth noticing that the neutralino density p^, evaluated according to Eq. (2), enters not 
only in C but also in ta (through C). Therefore the use of a correct value for p^ (rescaled 
according to Eq. (2), when necessary) is important also in determining whether or not the 
equilibrium is already set in a macroscopic body. 

Explicit calculations over the whole parameter space, given in the next sections, show 
that, whereas for the Earth the equilibrium condition depends sensitively on the values of 
the model parameters, in the case of the Sun equilibrium between capture and annihilation 
is typically reached for the whole range of m^, due to the much more efficient capture rate 
implied by the stronger gravitational field [25]. 

The annihilation rate given above refers to a macroscopic body as a whole. This is 
certainly enough for the Sun which appears to lis cLS db point source. On the contrary, in 
the case of the Earth, one has also to define an annihilation rate referred to a unit volume 
at point r from the Earth center: 
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TannI?^) = ^<av > n^{r) , (5) 
where n(r) is the neutrahno spatial density which may be written as [29] 

n(r) =noe-°"^^'''; (6) 
here a = 2nGp/ (3T) and no is such that 

Tann = ^<av> / d^rn\r) . (7) 



3.1 WIMP— nucleon cross section: hadronic uncertainties 

In Ref. [31] it is stressed that the couphngs between Higgs-bosons (or squarks) with nu- 
cleons, which typically play a crucial role in the evaluation of the neutralino-nucleus cross 
section, suffer of large uncertainties [32]. Actually, these couplings are conveniently ex- 
pressed in terms of three hadronic quantities, i.e.: the pion-nucleon sigma term 

o'nN = 2^^^ ~^ ^ N\uu + dd\N > , (8) 
the quantity do, related to the size of the SU(3) symmetry breaking, 

(To = ^{mu + rrid) < N\uu + dd- 2ss\N > , (9) 

and the mass ratio r = 2ms/ {i^u + ^d)- 

Because of a number of intrinsic theoretical and experimental problems, the determi- 
nation of these hadronic quantities is rather poor. Conservatively, their ranges can be 
summarized as follows (we refer to Refs. [31, 8] for details): 

41MeV<(T^;v <73MeV, (10) 



ao = 30 ^ 40 MeV , (11) 

and 

r = 29±7. (12) 

In the present paper, in order to display the influence of the uncertainties due to the 
hadronic quantities on the signals at neutrino telescopes, we will report our results for 
three different sets of values for the quantities (o",rAf , ctq ? '") as shown in Tab. 1. The set 
REF corresponds to the set of value referred to as reference point in Ref. [8] . The sets MIN 
and MAX listed in Tab. 1 bracket the range of hadronic uncertainties. 
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hadronic set 




ao[MeV] 


r 


MIN 


41 


40 


25 


REF 


45 


30 


29 


MAX 


73 


30 


25 



Table 1: Set of values for the hadronic quantities considered in the numerical analysis. 

In the case where the neutralino-nucleus interaction is dominated by the exchange of 
Higgs bosons, it is straightforward to estimate by how much the capture rate C is affected 
by the hadronic uncertainties. Indeed, in this case the dominant term in the interaction 
amplitude of the neutralino-nucleus scattering is provided by couphng between the two 
CP-even Higgs bosons and the down-type quarks: 

2 / 23 23 \ 

gd = ^ ifriN + —cr^N + —r {a^N - cro)j , (13) 

where m^r is the nucleon mass. Then: 

CmIn/C*REF — {gd,MlN / gd, REfY , CmAx/C*REF — igd,MAx/ gd,REFY ■ (14) 

Using the values of Table 1 for the three sets of hadronic quantities, one finds for g^. gd,MiN 
= 99 MeV, gd,REF = 290 MeV, gd,MAx = 598 MeV, respectively. We thus conclude that, 
because of the hadronic uncertainties, the capture rate in the case of set MIN is reduced 
by a factor ~ 9 as compared to the capture rate evaluated with the set REF, whereas C, 
evaluated with set MAX, is enhanced by a factor ~ 4. 

The consequences over the annihilation rate Fann is more involved, since the capture 
rate C enters in Fann not only linearly but also through r^. When in the celestial body 
capture and annihilation are in equilibrium (t > ta), one has Fann ~ C/2; then Fann, 
as a function of the hadronic quantities, rescales as C (see Eq. (14)); however, when the 
equilibrium is not realized, the uncertainties in Fann can be much more pronounced. For 
instance, for t <^ ta, Fann is proportional to C^, thus the rescaling factors for Fann are 
the squares of those in Eq. (14). 

These estimates will be confirmed by the numerical analysis displayed in the following 
section. 

3.2 Numerical evaluations 

The left panel of Fig. 1 shows the scatter plots for the ratios of the capture rates Cf^^^^/ C^ef^ 
(where i = set MIN, set MAX), obtained with the scanning of the parameter space illus- 
trated in Sect. 2. One sees that, as anticipated in the previous section, the numerical values 
accumulate (most significantly for light masses), around the numerical factors shown in 
Eq. (14). 
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The scatter plots for the ratios Tf^^^yr^^^^'^REp (where i = set MIN, set MAX) are 
displayed in the right panel of Fig. 1. As expected and discussed before, these numerical 
values are much larger as compared to those of the ratios Cf^'^^^/C^^p^, since many super- 
symmetric configurations are not able to provide a capture-annihilation equilibrium inside 
the Earth. This interpretation is validated by the scatter plots of Fig. 2, where one sees 
that for most configurations the equilibrium factor tanh^(to/'r4) deviate largely from the 
equilibrium value tanh^(to/'r4) = 1- 

The dependence of the annihilation rate for the Sun, on the hadronic uncertainties 
is shown in Fig. 3. Since the capture-annihilation equilibrium is realized in the Sun for 
all supersymmetric configurations of our model, one has here that r|j^j^ i/riNN ref — 
Cj^"°/C'ref) which implies that Ff^^ j^^x/rANN.REF ~ 4 and F^'g'j^ ^^jn/F|^n j^^p > 1/9. 
This is at variance with the case of the Earth which we have commented before. 

Moreover, one notices from Fig. 3 that for many supersymmetric configurations F^^j^ 
depends very slightly (or negligibly) on the variations in the the hadronic quantities. This 
is due to the fact that on many instances the capture of neutralinos from the Sun is 
dominated by spin-dependent cross-sections. 



4 Calculation of neutrino fluxes 
4.1 Neutrino production 

Once the Dark Matter particles are accumulated in the center of the Sun and the Earth, 
they can annihilate into leptons, quarks, gauge and higgs bosons. These particles then 
decay producing neutrinos. To calculate the neutrino spectra coming from each annihi- 
lation channel, lepton decay and quark hadronization process have to be considered. In 
[33], the authors have used a PYTHIA Monte Carlo simulation to calculate the spectra 
of neutrinos, coming from Dark Matter annihilation in the Sun and in the Earth, for the 
following channels: bb, rf, cc, qq, gg (with q = u, d, s quarks). Three main differences and 
improvements have been implemented in [33] with respect to previous calculations. The 
first one is the prediction of the neutrino spectra for the different neutrino flavors: z/g, 
and i>r (not only i/^, like in previous works). The second main improvement consists in an 
appropriate implementation of the energy loss that hadrons and leptons can undergo be- 
fore decaying. Finally, the third difference is represented by the calculation of the neutrino 
spectra for light quarks u, d, s, that were usually neglected in previous calculations. 

We use the initial neutrino spectra reported in [33] for bb, tt, cc, qq and gg annihilation 
channels. The annihilation of two neutralinos can also produce two higgs bosons or one 
gauge and one higgs boson in the final state, although these two channels (as well as the 
annihilation channels into tt and into two gauge bosons) are absent for < 80 GeV (we 
remind that our model has an absolute lower limit on the higgs mass of 90 GeV). 
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4.2 Neutrino propagation 

For a precise estimate of the neutrino flux at tlie detector site, it is important to take into 
account the main processes that can occur during the neutrino propagation: the oscillation 
and the incoherent interaction with matter. These effects have been vastly analyzed in 
[33, 34], and then applied to specific model-dependent studies, see e.g. [35]. 

The equations that describe the evolution of the neutrino spectra can be formally 
written using the density matrix formalism: 



_l_ dp 

NC d^ 



(15) 

cc 



The first term describes the oscillation of neutrinos in matter, where the Hamiltonian is 
given by the sum of the mass matrix in the weak basis (z/g, z/^, v-r) and the Wolfenstein 
potential: 

H^ = ^±V2GfN, diag(l, 0, 0) (16) 

holds for (anti-)neutrinos for the minus (plus) sign. We fix the values for the neutrino 
mixing angles and squared mass differences to the best fit values reported in [36], and 
we set the mixing angle ^13 to zero. A different choice of ^13 would marginally affect the 
prediction on the neutrino flux, as reported in [33, 34]. 

The second term in Eq. (15) takes into account the neutrino energy loss and the rein- 
jection due to neutral current interactions. The last term, instead, represents the neutrino 
absorption and Vr regeneration through charge current interactions. For the explicit form 
of each term and for exhaustive explanations we refer to [33]. 

For Dark Matter annihilation inside the Sun, the integro-differential equation for the 
density matrix has to be solved numerically to find the neutrino spectra at the detector 
site. For simplicity, we neglect the z/^ regeneration effect, since it provides only a negligible 
correction for the WIMPs mass range of our interest: m-^ < 80 GeV. 

In the case of annihilations inside the Earth's core, the calculation of the neutrino 
spectra can be further simplified. Indeed, the interactions with matter can be neglected, 
since the mean free paths of neutrinos, in the core and in the mantle, are much bigger than 
the Earth's radius i?® for < 10 TeV (for anti-neutrinos the mean free paths are almost 
a factor two greater, due to the difference in the cross-sections): 

y^ore ^ _^ _ g^g ^ 10 \ „ ^® . , A"^""*'^ = ^} „ ~ 9.2 X 10^ 



a^N^""-^ ■ (E^/GeV) ' ^^Nmantie ■ (E^/GeV) ' 

(17) 

Therefore, for the propagation inside the Earth, it can be safely taken into account only the 
oscillation effect. Moreover, for E^ > 1 GeV, the dependence on the "solar" parameters 
Amgi and 612 is extremely weak and can be neglected. Since we are considering vanishing 
^13, Earth's matter effects are negligible and the neutrino oscillation is driven by the 
"atmospheric" parameters Am|^ and 623. In this case, the main oscillation channel is the 
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i^fi ^ i^T and the value of the oscillation and the survival probability Pafs is simply given 
by the vacuum two-flavors formula: 



P^p{r, E,) = S^fs - e^p sm\2e) sin^ M .27 ^ /GeV) ) 



where the parameter is equal to 1 (-1) for a = /3 (a 7^ 

In the case of the Earth, the differential muon-neutrino flux at the detector, as a function 
of the zenith angle 6'^, can be written as: 

where the function GapiOz, Ey) encodes the dependence on the oscillation probability and 
on the WIMP distribution inside the Earth. Using Eq. (6), we find the following expression: 

G^p{9z, E,) = ^ / dr exp [-2 a (r^ + i?^ - ry)] P^p{r, E,) , (20) 

with y = 2 i?® cos6'„. On = re — 6z and /3 = a i?^. The differential muon anti-neutrino flux 
at the detector can be obtained by a formula analogous to Eq. (19). 



5 Muon fluxes 

For the calculation of up-going muons we follow the formalism described in [37], to which 
we refer for specific details. In these references, it has been shown that the differential 
muon flux is given by the following expression: 

= -J— r,E. r '-^^^^^ . (21) 

d COS 6 zdEf, a + bEf,JE^ J ^ d cos 9 ^dE^ dE'^ 

The quantities a and h parametrize the energy loss due to ionization and to radiative effects, 
Na is the Avogadro's number and da^^/dE'^ is the differential charge current cross section, 
which is mainly due to deep inelastic scattering, for energies E^ > 1 GeV. The total muon 
flux is then divided in through-going and stopping muons: 

*f (cos^^) = r dE,-^^A'^^{L{E,\ez). (22) 

^ A{Lmin,Oz) Jeii^ dcosOzdE^ 

with Ejl^ being the energy threshold of the detector for upward-going muons, L the muon 
range in water ( L = ^ In '^^l^^ ) and Lmin = L{E^^). The function A^''^{L, 9z) represents 
the effective area for through-going and stopping muons, while, A{Lminidz) is the total 
effective area of the detector, i.e. the projected area that corresponds to internal path- 
lengths longer than Lmin-, for a fixed value of the zenith angle 9z- 
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The classification of upward-going muons into tlie two subcategories reported above is 
strictly detector-dependent, since it depends on the shape and the size of the detector. For 
a detector with cylindrical geometry (with radius R and height H), it has been shown in 
[38] that the function A{L, 9z) acquires the form: 



= 2RHsixi9^^l-x^ + 2i?2|cos^^| cos"^ x - 3a;Vl - x"^ Q{Lmax{Oz) - L) , 

(23) 

with X = Lsm6z/2R and Lmax{6z) = min [2i?/ sinOz, H/\cos6z\]. The effective area for 
stopping muons A^{L, 9z) is given by Eq. (23) and the one for through-going is A'^{L, 6^) = 

[A{Lrmn,ez)-Amez)]. 

We will focus our analysis to the SK detector. Indeed, even if the recent results of 
Ice-Cube 22-strings [39] improve the SK bound in the high mass region (m^ > 200 GeV 
for the hard channel and > 500 GeV for the soft channel), they do not constrain 
the parameter space at low masses. The Ice-Cube detector, provided with the Deep-Core 
arrays [40], will improve significantilly the SK bound for mass > 40 GeV, as has been 
shown in [41], but not for the low mass range in which we are mainly interested: 7 GeV 
< < 80 GeV. 

The SK detector has R = 16.9 m, H = 36.2 m and an energy threshold Ejj^ = 1.6 GeV 
(that corresponds to Lmin = 7 m). Using these features the effective area can be calculated 
through Eq. (23) and, with Eq. (22), the muons coming from Dark Matter annihilation can 
be divided in through-going and stopping muons. The same formalism can be applied to 
the calculation of the expected muons background coming from atmospheric neutrinos. 
We use the atmospheric neutrino flux from Honda et al. [42] and we reproduce with great 
accuracy the zenith angle distribution for stopping and through-going muons, as predicted 
by the SK collaboration. For the muon energy loss in rock and in water, we use the 
tabulated values reported in [13]. Following the analysis for through-going muons of [14], 
the SK limit on the muon flux can be defined as: 

m; 90% C.L.) = (24) 

where Ni^q is the upper Poissonian limit at the 90% C.L., given the measured events and 
the muons background from atmospheric neutrino, and T is the detector lifetime. We 
do not consider the detector efficiency, since it is almost equal to 100% for upward-going 
muons. Using Eq. (24) and the SK data collected from May 1996 to July 2001 [45], we find 
the following limits on through-going ($^) and stopping ($^) muons: 

^l,Earth ^ 0.8 X 10"^^ cm"^ s~^ at 90%C.L. , (25) 

^lEarth < 0.5 X 10-14 cm-2 at 90%C.L. , (26) 
"^Isun ^ 1-2 X 10-^4 cm-2 s-^ at 90% C.L. , (27) 
K,sun < 0.5 X 10-14 cm-^s-i at 90% C.L. , (28) 
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The limits for the Earth are obtained considering the angular bin —1.0 < cos^^ < —0.9, 
while the values reported for the Sun are the average limits obtained varying the zenith 
angle from cos9z = —1.0 to cos9z = 0. 

The values reported in Eqs. (25)-^(28) have to be compared with the muon flux induced 
by neutralinos annihilation. For the case of the Earth, we fix the angular opening to 
— 1.0 < cos 6*2 < —0.9 also for the calculation of muons coming from Dark Matter. For the 
case of the Sun, we divide the muons in stopping and through-going using the SK effective 
area averaged over the zenith angle. 

In the calculation of the muon flux, we neglect the kinematical angle between the 
neutrino and muon direction, which can be relatively large for muons close to threshold. 
In any case, the average deflection angle is at most of the same order of the angular bin 
over which we integrate our signal, for the stopping and through going muons. Considering 
also the detector resolution, to neglect the kinematical angle does not affect our results in 
a relevant way. This is conflrmed by the quite good agreement we obtain in our calculation 
of the atmospheric-neutrino events with the SK evaluation [46]. 

6 Results and conclusions 

In this section we give our results for the muon fluxes expected at a neutrino telescope 
with a threshold muon energy of 1.6 GeV, generated by annihilation of light neutralino 
pair-annihilation in the Earth and in the Sun. 

6.1 Fluxes from the Earth 

The upper panel of Fig. 4 displays the scatter plots for the expected muon flux integrated 
over the muon energy for i?^ > 1.6 GeV for the upward through-going muons. The three 
columns refer to the evaluation of the fluxes using in turn the three different set of hadronic 
quantities deflned in Sect. 3.1. 

The main features common to the three scatter plots are easily interpretable in terms 
of the following properties: 

i) The various peaks for < 40 GeV are due to resonant capture of neutralinos on 
oxygen, silicon and magnesium; indeed, these elements are almost as abundant in Earth 
as iron, which is the most relevant target nucleus for the capture of neutralinos of higher 
mass. The dip at ~ 45 GeV is a consequence of a depletion of the neutralino local 
density, implied by the rescaling recipe of Eq. (2) and a resonant effect in the (Z-exchange) 
neutralino pair annihilation when < mz/^ (note that the neutralino relic abundance 
is inversely proportional to the neutralino pair-annihilation). 

ii) Also the fact that the muon signal for light neutralinos {m^ < 25-30 GeV) is lower 
than the one at higher masses can be understood. Indeed, for light neutralino masses the 
branching ratio of the annihilation process into the r — f flnal state, which is the one with 
the highest neutrino yield per annihilation, is suppressed. This last property being in turn 
due to the fact that, for these masses, the flnal state in 6 — 6 in the annihilation cross 
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section has to be the dominant one in order to keep the neutrahno rehc abundance below 
its cosmo logical upper bound [1]. 

iii) Moreover, lower masses imply softer neutrino spectra which entail fewer muons 
above threshold. 

The comparison of the fluxes in the three columns shows how relevant can be the role 
of the size of the hadronic quantities on the final outputs. The suppression (enhancement) 
of the flux in the case of the set MIN (MAX) as compared to the flux for the set REF are 
set by the numerical factors previously discussed for T^^^. This entails that, whereas the 
overall muon flux is completely below the present experimental bound in the case of the 
minimal set of the hadronic quantities, some part of the spectrum would emerge sizably 
above the limit for neutralino masses > 50 GeV for the other sets. In the case of set 
MAX, neutralino configurations with masses ~ 15 GeV or ~ 25-30 GeV might 
produce some measurable signal. 

Since also the dependence of the muon signals on the astrophysical parameters Vq and 
Po is important, in the lower panel of Fig. 4 we display the through-going fluxes for the 
three representative values of Vq and po which we discussed in Sect. 3. The overall increase 
in the fluxes in moving from left to right is essentially due to the increase in the value of 
the local DM density. In these scatter plots the hadronic quantities are set to the value 
REF. 

The fluxes for upward stopping-muons from the Earth are given in Fig. 5. The scheme 
of this figure is the same as the one of the previous Fig. 4: the dependence of the fluxes 
on the hadronic quantities can be read in the upper panel, the one on the astrophysical 
parameters is displayed in the lower panel. 

Because of the uncertainties affecting the evaluations of the muon fluxes, mainly due 
to the hadronic quantities, we cannot convert these results in terms of absolute constraints 
on supersymmetric configurations. However, we can conclude that the analysis of stopping 
muons from the Earth can have an interesting discovery potential not only for masses 
above 50 GeV, but also for light neutralinos with ~ 15 GeV or ~ 25-30 GeV. 
Notice however that the neutralino configurations which provide the highest values of the 
muon fluxes, mainly at ~ 50-70 GeV, are actually disfavored by measurements of 
WIMP direct detection [ 17] which have their maximal sensitivity in this mass range. 

6.2 Fluxes from the Sun 

The fluxes of upward through-going muons and for stopping muons from the Sun are 
provided in Fig. 6 and in Fig. 7, respectively. The schemes of these figures is the same as 
the ones of the previous Fig. 4 and Fig. 5, respectively: the dependence of the fluxes on the 
hadronic quantities can be read in the upper panels of Fig. 6 and Fig. 7, the one on the 
astrophysical parameters is displayed in the lower panels of the same figures. 

From these results one notices that through-going muons can only be relevant for neu- 
tralinos with masses > 50 GeV or ~ 35-40 GeV, whereas stopping muons can 
potentially provide information also on some supersymmetric configurations with masses 
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down to ~ 7 GeV, in the favorable cases of high values of the hadronic quantities and 
of the astrophysical parameters. 

6.3 Fluxes of stopping muons for supersymmetric configurations 
selected by the DAMA/LIBRA annual modulation data 

Now we give the expected upward muon fluxes from Earth and Sun which would be pro- 
duced by neutralino configurations which fit the annual modulation data of the DAMA/LIBRA 
experiment [6]. As before, for definiteness the analysis is performed in the framework of 
the isothermal sphere. The selection of the supersymmetric configurations is performed on 
the basis of the analysis carried out in Ref. [8]: for any set of astrophysical parameters and 
hadronic quantities, from the whole neutralino population are extracted the configurations 
which fit the experimental annual modulation data, and the relevant muon fluxes are eval- 
uated. As for the yearly modulation data, we consider both outputs of the experimental 
analysis of the DAMA Collaboration: those where the channelling effect [5] is included as 
well as those where this effect is neglected. We recall that the way by which the channeling 
effect has to be taken into account in the analysis is still under study; thus the actual phys- 
ical outputs in the analysis of the experimental data in terms of specific DM candidates 
could stay mid-way, between the case defined as channeling and the no-channeling one, 
respectively. 

We only report the results for stopping muons, since, as we have seen above, this is 
the category of events which can provide the most sizable signals. To avoid proliferation 
of figures, only fluxes calculated with the set REF for the hadronic quantities is reported 
here. 

Fig. 8 displays the fluxes for the upward stopping muons expected from the Earth in 
case of no-channeling (upper panel) and in the case of channeling (lower panel). The 
corresponding fluxes from the Sun are shown in Fig. 9. 

We note that depending on the role of channeling in the extraction of the physical 
supersymmetric configurations, the stopping muon fluxes can have a discovery potential 
with an interesting complementarity between the signals from the two celestial bodies: 
whereas the flux from the Earth cannot give insights into neutralino masses below about 
15 GeV, the flux from the Sun would potentially be able to measure effects down to ~ 
7 GeV. 

It is worth remarking that under favorable conditions provided by the actual values of 
the involved parameters, a combination of the annual modulation data and of measure- 
ments at neutrino telescopes could help in pinning down the features of the DM particle 
and in restraining the ranges of the many quantities (of astrophysical and particle-physics 
origins) which enter in the evaluations and still suffer from large uncertainties. 

We stress once more that the present analysis, for definiteness, was performed only in 
the standard case of a halo DM distribution function given by an isothermal sphere. Use 
of different halo distributions such as those mentioned in Sect. 3 could modify sizably the 
role of specific supersymmetric configurations. 
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We wish here to recall that indirect signals of light neutralinos could also be provided 
by future measurements of cosmic antideuterons in space [10]. Finally, investigations at 
the Large Hadron Collider will hopefully provide a crucial test bench for the very existence 
of these light supersymmetric stable particles [ lH] . 
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Figure 1: Ratios of capture rates (left panel) and annihilation rates (right panel), in the 
case of the Earth, calculated for the hadronic sets MIN and MAX with respect to the 
hadronic set REF. The local rotational velocity is set to its central value: Vq = 220 Km 
s-i (po = 0.34 GeV cm'^). 




Figure 2: Equilibrium factor tanh^ (Iq / ta) , in the case of the Earth, displayed for the 
hadronic sets MIN, REF and MAX. The local rotational velocity is set to its central value: 
vo = 220 Km s^^ (po = 0.34 GeV cm'^). 
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Figure 3: Ratios of annihilation rates, in the case of the Sun, calculated for the hadronic 
sets MIN and MAX with respect to the hadronic set REF. The local rotational velocity is 
set to its central value: = 220 Km s~^ (po = 0.34 GeV cm~'^). 
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Figure 4: Upward through-going muon flux, generated by hght neutrahno pair-annihilation 
inside the Earth. The upper panel shows the dependence of the muon flux on the hadronic 
quantities, for fixed values of the astrophysical parameters: Vq = 220 Km s~^ and po = 0.34 
GeV cm""^. The lower panel shows the dependence of the muon flux on the local rotational 
velocity Vq and the total DM density po, for the hadronic set REF. The horizontal line 
represents the experimental limit on through-going muons from the Earth obtained using 
the SK data, see Eq. (25). 
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Figure 5: The same as Fig. 4, but in the case of upward stopping muons. In this case, the 
horizontal hue refers to the experimental limit on stopping muons from the Earth obtained 
using the SK data, see Eq. (26). 
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Figure 6: The same as Fig. 4, but in the case of hght neutrahno pair-annihilation inside 
the Sun. In this case, the horizontal line refers to the experimental limit on through-going 
muons from the Sun obtained using the SK data, see Eq. (27). 
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Figure 7: The same as Fig. 4, but in the case of hght neutrahno pair-annihilation inside 
the Sun and of upward stopping muons. In this case, the horizontal line refers to the expe- 
rimental limit on stopping muons from the Sun obtained using the SK data, see Eq. (28). 
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Figure 8: Upward stopping muon flux, generated by light neutralino pair-annihilation 
inside the Earth. The configurations displayed are only the ones compatible with the 
DAMA/LIBRA annual modulation region, obtained without including the channeling effect 
(upper panel) and including the channeling effect (lower panel). The three columns show 
the results for the different sets of astrophysical parameters, defined in Sect. 3. The set 
REF is used for the hadronic quantities. The horizontal line represents the experimental 
limit on stopping muons from the Earth obtained using the SK data, see Eq. (26). 
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Figure 9: The same as Fig. 8, but in the case of hght neutrahno pair-annihilation inside the 
Sun. In this case, the horizontal line refers to the experimental limit on stopping muons 
from the Sun obtained using the SK data, see Eq. (28). 
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